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We assess the relative contribution to the COBE - measured microwave anisotropy arising both 
from relic gravity waves as well as primordial density perturbations originating during inflation. We 
CO 1 show that the gravity wave contribution to the CMBR anisotropy depends sensitively upon n - the 
primordial spectral index (\5k\ 2 oc fc n ), increasing as n deviates from a Harrison - Zeldovich spectrum 
(n = 1). As a result, for n < 0.84 the contribution from gravity waves towards ST/T is greater than 
the corresponding contribution from density perturbations, whereas for n > 0.84 the reverse is 



i 



Q_i[ true, (n — 0.84 corresponds to an expansion index p = 13.5 in models with power-law inflation 
a ex t p . ) Our results show that for a scale- invariant Harrison-Zeldovich spectrum generated by 
• • chaotic inflation, gravity waves contribute approximately 24% to the CMBR anisotropy measured 
. ^ by COBE. Applying our results to the cold dark matter scenario for galaxy formation, we find 
that in general CDM models with tilted power spectra (n < 1), require the biasing parameter 
^ to be greater than unity, on scales of 16/i^g 1 Mpc. We also obtain an expression for the COBE - 
normalised amplitude and spectrum of the stochastic gravity wave background and compare it with 
the sensitivity of planned laser-interferometer gravity wave detectors. 
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Both gravity waves and adiabatic density perturbations, characterised by a power spectrum 
\Sk\ oc k n , are generically predicted by the Inflationary Universe scenario. On entering our horizon 
during the epoch of matter-radiation decoupling, density perturbations as well as gravity waves cause 
distortions in the cosmic microwave background radiation (CMBR) due to the well known Sachs- 
Wolfe effect. The anisotropy in the CMBR recently observed by COBE measures the amplitude of 
gravity waves as well as density perturbations providing thereby a direct probe of the inflationary 
epoch of our universe. 

The close similarity between gravity waves and massless scalars, in a Friedman - Robertson 
- Walker (FRW) background originally pointed out by Grishchuk 1 , makes it possible to express 
each polarisation state of the graviton in terms of solutions to the massless, minimally coupled 
Klein-Gordon equation: /i x (+) = V&ttG 0^ e _lk x 



Xk + [k 2 - -] Xk = (la) 

where \k = <Pk x a(v)> a(v) is the scale factor of the Universe: ds 2 = a 2 (rj)(di] 2 — dx 2 ), and k 
is the comoving wavenumber k = 2ira/\. Eq.(la) closely resembles the Schrocdinger equation in 
quantum mechanics, with a/a, playing the role of the potential barrier 'V. (The form of V = (a/a), 
is shown in Figure 1, for inflation followed by a matter dominated epoch.) Solutions of equation (la) 
have some interesting properties, for instance, small wavelength solutions of (la) are adiabatically 
damped: 

4>t{r}) = -^r = r-^ - exp(-ifc?7), (16) 

rfcw/ ^1]) fcr,»27rV2fca(?7) 

whereas long wavelength modes asymptotically approach the form : 

Mv)^^ = A(k)+B(k) , (lc) 

a(rj) fcr,«27r J a 2 

where A and B are constants. The form of (lc) implies that the amplitude of a fluctuation (equiva- 
lently - gravity wave), freezes to a constant value when its wavelength (during inflation), becomes 
larger than the corresponding Hubble radius. Since modes are continuously being pushed outside 
the Hubble radius during inflation, a mode with a fixed wavenumber "k", which left the Hubble 
radius during inflation (k = 2ir r/f : ), will re-enter it later during the epoch of matter domination 
(k = 2n tjq 1 ). The amplitude of the mode between these two epochs is supcr-adiabatically amplified, 
which in field-theoretic language corresponds to the quantum creation of gravitons 1 (see fig. 1). 

The possibility that relic gravity waves may be generated during inflation, was first investigated 
by Starobinsky 2 , who showed that the spectrum of gravitons is scale invariant in a fairly broad 
wavelength interval, if inflation proceeds exponentially and is followed by a radiation dominated 
epoch. The quantum creation of relic gravity waves was subsequently extended to other inflationary 
models by Abbott & Wise 4 and Sahni 6 . (The effect of a post inflationary matter dominated epoch on 
the graviton spectrum was studied by Allen 5 , and also by Sahni 6 .) In particular it was demonstrated 
that, for power law expansion a = (t/t ) p = ("^)^ — i { u ~ § = j^j = — , where w = ^ is the 



equation of state of matter) Eq.(la) can be solved exactly, with the following solution describing 
positive frequency waves in an FRW Universe in the adiabatic limit: 



(y = I corresponds to exponential inflation, w — — 1; ^ > | to power law inflation, — 1 < w < 
— i;i<^<|to super-exponential inflation (pole - driven inflation), w < — 1; and — | < v < to 
matter dominated expansion, < w < 1 .) From (lc) and (Id) it can be shown that, the amplitude 
of gravity waves on scales greater than the Hubble radius, has the time independent form 6 



V27T fe| m |=27T \ m pJ 



A{v)=Y( V )^{v-\r\ 



(2) 



where we have summed over both polarisation states of the graviton. The value of u is related to the 
spectral index n of density fluctuations as v = (4 — n)/2 and Hi is the value of the Hubble parameter 
at a time k\r]i\ — 2tt, when the given mode left the Hubble radius. For exponential inflation v = | 
and h+ tX = (^/V^)(Hi/m p ). (m p is the planck mass, m p = 1.2 x 10 19 GeV.) 

Both long wavelength gravity waves as well as adiabatic density perturbations generated dur- 
ing inflation, create potential fluctuations at the surface of last scattering, when their wavelength 
becomes smaller than the horizon size during the matter dominated epoch. This leads to distor- 
tions in the cosmic microwave background, described by the well known Sachs - Wolfe effect 4,7 
(ST/T) ~ | (Sf/c 2 ). Writing (ST/T) in terms of a multipole expansion: 

^ = ]Ta im YT(M), (3) 



T 

l.m 



it can be shown that, for gravity waves 



4.8 



a| = (|a 2m | 2 ) =0.145 h 2 +>x g(n), g(n) = 5770.3 x jf db ( -J [7 2 (6)] 2 , (4a) 



where 8 



' lh: L, . — ^r^r- b=k ^- {Ah) 

We have numerically integrated (4b) to obtain g(n) {kr\ rec w 0), which is well approximated by the 
fitting formula g(n) = cxp[0.6(n — 1)] over the range 0.5 < n < 1. The corresponding value of the 
rms quadrupole amplitude is 

Q 2 T = i L 4F 2 = 2.9 A 2 {u) (^-) g{n) F 2 . (4c) 



4-7T \ m \ , 

(For exponential inflation v = |, and = (2.9/47T 2 )(-Hi/m p ) 2 F 2 .) (The subscript "T" , denotes the 
fact that the quadrupole anisotropy is induced by tensor waves, gravity waves being quadrupolar, do 
not contribute to a dipole component ai.) F 2 incorporates the finite beam width of the COBE-DMR 
instrument (F 2 f=s 0.99) . 

Scalar density perturbations generated during inflation, are characterised by a spectrum \Sk\ 2 = 
A k n , the scale invariant Harrison - Zcldovich spectrum n = 1, is predicted by inflationary models 
with exponential expansion. Models with power law inflation (a(t) oc t p ,p > 1), which arise in 
a number of theories including extended inflation 9 , predict a more general value for the spectral 
index 10 , n = (p — 3)/(p — 1) < 1 . (n > 1 , is predicted for super-exponential inflation 28 .) 



The amplitude of density perturbations at horizon crossing (k — Hq/c) is given by 11 



p J H 2ir 2 2tt 2 

(we assume c = 1, for simplicity.) The associated fluctuations in the microwave background are 
given by ( for I > 2) 

aN(|a im | 2 > = g j ( 00 §|4| 2 . ? f(M (6a) 

where x = (2/H ), is the present day horizon size (jf(kx) are spherical Bessel functions). For power 
law spectra |<5fc| 2 = A k n , we obtain 12 



AH^ ^ fSp\ 2 r(3-n)r(3±a) 



Values of n in the range 0.5 < n < 1.7, are consistent with the recent COBE results 13 . For n = 1, 
(6b) reduces to a 2 . = (tt/12) (8p/p) 2 H . 

The related value of the rms quadrupole amplitude is given by 

(the subscript "S" refers to scalar perturbations.) 

As shown by a number of authors, scalar field fluctuations which left the Hubble radius during 
inflation, upon re-entering the horizon during radiation (matter) domination, give rise to a density 
fluctuation, whose rms value is given by 14 

'A =^ (8) 

where Hi is the Hubble parameter at a time t\ 1 when a scale which reenters the horizon at t , left the 
Hubble radius during the inflationary epoch (see Figure 1). For modes entering the horizon during 
radiation (matter) domination, b = 4(|). (The main contribution to the observed quadrupole 
anisotropy however, comes from modes entering the horizon during matter domination, we shall 
therefore assume b — | in the ensuing discussion.) 6<fr arises because of quantum fluctuations in the 
scalar field during inflation. For modes entering the horizon today 5(j) is given by 5<fi — A(v)H\ (see 
(2)), from which we recover the standard result 5(f) — for exponential inflation 14 ' 15 . 

For power law inflation, the inflaton potential has the form V(<p) = 
V exp(— y/Wir/p (0/m p )) which results in the following exact solution 16 to the field equations: 
H = (p/t), (H/<j)) = (y/4Trp/m p ), as a result 



Combining Eqs. (7) and (9), the rms quadrupolc amplitude of the temperature anisotropy 



H^ 2 



Q*=- p A*{v)[-± \ f(n)F 2 (10) 
10 \m p J 

(y — | = (p — 1) _1 = (1 — n) /2) ), for power law inflation. (See also Fabbri, Lucchin and Mataresse; 
Lyth & Stewart.) 17 The above result for power law inflation can be extended to cover models of slow- 
roll (quasi-exponential) inflation . As shown in Souradeep and Sahni 18 , Q 2 S for a chaotic inflation 
model with a m 2 (j) 2 potential corresponds to a power law model with p — 2N, where ./V = H At is the 
number of e-foldings (At is the duration of the inflationary epoch and H w constant is the Hubble 
parameter during inflation) . In terms of the number of c- foldings in a model of slow-roll inflation, 
the following relation 

2N = p=^, (11) 
n — 1 

fixes the spectral index of density perturbations n. As a result slow-roll inflation ( chaotic m 2 <p 2 
model) with minimal e-foldings N — 67 is equivalent to a power law model with p = 134 which 
corresponds to a spectral index n — 0.985. 

Since the coefficients &i m in (3) are Gaussian random variables, the variances Q\ and Q 2 S given 
by (4c) and (10) must be added, when considering the combined effect of both gravity waves and 
adiabatic density perturbations on the anisotropy of the CMBR. As a result, the net quadrupolc 
anisotropy will be given by 

Qcobe-dmr = Qs + Qt (12) 

where Qcobe-dmr. is the rms-quadrupole normalised amplitude Q rm s-ps = 5.86 x 10^ 6 detected 
by COBE 13 . From (4c) and (10), we find that the ratio (Q%/Q%) = (2.94 (n - l)g(n)/(n - 3)/(n)) 
depends sensitively on n the spectral index of density perturbations. We find that for n > 0.84, the 
contribution from gravity waves to the CMBR anisotropy dominates the contribution from scalar 
density perturbations, whereas for n < 0.84 the reverse is true, (n — 0.84 corresponds to p = 13.5.) 
The ratio (Qx/Q's) nas been plotted against n in Figure 2 . For slow-roll (chaotic) inflation we find 
that the microwave distortion due to gravity waves amounts to ~ 31% of the distortion caused by 
density fluctuations, in agreement with earlier work by Starobinsky 8 and contrary to recent claims 
by Krauss and White 19 . 

Substituting the values of Q\ and Q 2 S in Eq.(12), we obtain 



- r 7r 2 n-3 



2 QcOBE-DMR 

10 n-l 



f(n) + 2.9 g{n) 



[A{v)]~ 2 (13) 



for power law inflation. Setting p = 2N in Eq. (13) we can recover the result for slow-roll inflation 18 . 
As a result Eq.(13) allows us to determine the Hubble parameter during inflation both for power 
law, as well as for quasi-exponential inflation. 

The value of (Hi/m p ) allows us to determine the value of Q 2 S - the fraction of the quadrupole 
anisotropy contributed by density perturbations, as well as (dp/ p) 2 H - the amplitude of density 
fluctuations at horizon crossing (see (5), (9)). The value of (Sp/p)n fixes the normalisation A of the 
power spectrum of density fluctuations to be 



A 

2tt 2 



72.5 n - 1 , x , 5tt 
l&T ^3 + 32 /W 



if, 



n+3 n 2 



COBE-DMR 



(14) 



( we have reintroduced c in the above expression for completeness and as noted earlier (Eq. 4c), 
the finite beam width factor F 2 ~ 0.99 for COBE-DMR ). We now use the COBE-DMR results to 
normalise the spectrum and predict the value of the rms mass fluctuation on a given scale 



M 



) {R) = i/ k2dk \ S k\ 2 W 2 (kR), 



(15) 



where W(kR) is the top hat window function 



20 



Our results for cold dark matter models with power law primordial spectra : \Sk | 2 = A k n Tfe iC( jm, 
(where Tk^dm is the CDM transfer function given in Appendix G of Bardeen et al. 21 ) are shown in 



Figure 3, for different values of the parameters n and (f! c dm — 0.9, ribaryonic — 0.1, is assumed). 
(Values of n in the range 0.5 < n < 1.7, are consistent with the COBE - DMR results 13 .) For a 
scale- invariant spectrum n = 1, our results agree with those of Bond and Efstathiou 22 . Our results 
also show that ^f(16h§ Mpc) < 1 for n < 0.87. Since mass need not trace light, in models with 
nonbaryonic dark matter, our results can be taken to mean that a biasing factor (&i 6 = 
mostly greater than unity, is required in order to reconcile theoretical models based on power law 
inflation, with observations. From Figure 3 it follows that the value of &16 is sensitive to both n 
as well as /150, decreasing with increasing n and /150. It may be noted that since gravity waves 
contribute predominantly to the CMBR anisotropy for n < 0.84, the biasing factor b 16h -i is much 
larger than it would have been had only the contribution from scalar density perturbations to the 
CMBR been considered. This makes CDM models with power law inflation less compatible with 
the excess galaxy clustering observed in the APM survey, than had previously been assumed 23 . 

Having normalised the amplitude of scalar density perturbations generated during inflation, we 
now proceed to evaluate the normalised spectrum of relic gravity waves. During inflation (ie for 
77 < 770 < 0), gravitons are described by Eq.(ld) which represents the state corresponding to the 
adiabatic vacuum . During the matter dominated epoch however, gravitons are described by a state 
which is a linear superposition of positive and negative frequency solutions of Eq.(la): 



Mv) = + 04t ] (16a) 

where 

r] > \t]q\. /i can be related to the equation of state of matter fj, = | j^^; where w = P/e, (— | < /j, < 
for reasonable equations of state for matter: < w < 1). 

The Bogoliubov coefficients a and (3 can be obtained by matching (Id) and (16) at wavelengths 
larger than the horizon size via (lc). As a result we obtain for krj < 2tt, 



a±(3 = ±i 1 ±1 (*^)=F("+H) (17) 

where 7 = 7 r- 1 r(^)r(l + and |a| 2 - |/?| 2 = 1, (sec Salmi 6 (1990) for details). (For kr^ > 2tt 
the adiabatic theorem gives a ~ 1, (3 ~ 0.) 

The energy density of created gravitons can now be determined exactly by substituting the 
expression for (j>k(v) m (16), with the values of the Bogoliubov coefficients given in Eq.(17) in the 



expression 29 



( T o) = ^2 I dkk "[ 4>k +fc 2 !^| 2 ), (18) 



27r 2 a : 

(both polarisation states of the graviton have been included in Eq.(18).) 



Of greater relevence to us is the spectral energy density of gravity waves, e(u>) — which 



can be derived from Eq.(18) and has the simple form 6 



e(w) = b\ u 1 - 2 *( — ) 2 e m for 1< Q < Afi; 



47T 



l2 -3-25 1 Hi ,2, f . ^ 



e(w) = 6fo) d - 2l/ ( — ) 2 e r for — tl r 3 < Co (19) 
m p 47T 

where Co = = is the dimensionless wavenumber expressed in units of the horizon scale (A 
being the physical wavelength and \h being the present scale of the horizon, A/j ~ 2. lO 28 /^ 1 cm ). 
I < Co < ^fl r 2 corresponds to wavelengths larger than the horizon size at matter radiation equality, 
e m and e r are the energy densities of matter and radiation respectively; e m w e cr = 4.2.10 _9 /i 2 
ergs/cm 3 , f2 r = ^- ~ 10 _4 /i^" 2 , /150 is the present value of the Hubble parameter in units of 50 
km/sec/Mpc. 

From Eq.(19) we see that the spectrum of the gravity wave background is uniquely specified 
once (— ) - the dimensionless value of the Hubble parameter during inflation is established from 
Eq.(13). The ratio of the spectral energy density of gravity waves to the critical energy density 
51 g (A) = e(X)/e cr and the dimensionless amplitude of gravity waves, /i x+ (A) have been plotted 
in Figure 4, for gravity waves created during exponential inflation, slow-roll (quasi-exponential) 
inflation and power law inflation. (For the corresponding analytical expressions describing O g (A), 
see Sahni 6 and Grishchuk & Solokhin 24 .) The quantities h x + (X) and £1 9 (A) are directly related 
by ft. 2 + (A) = 24(A/A/ l )£l 9 (A). We find that the amplitude of gravity waves is significantly smaller 
than the sensitivity of the current generation of terrestrial bar and beam detectors 25 . The best 



hope for the detection of the stochastic gravity wave background appears to lie with the space - 
interferometer, whose development is still in the conceptual stage 26 . 

After completing this work, we became aware of papers 27 by Davis et al.; Salopek; Liddle & 
Lyth; Lucchin, Mataresse & Mollerach ; and Lidsey and Coles, in which results overlapping with 
those of this paper and ref [18] were obtained. 
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Figure Captions 



Fig. 1 

The superadiabatic amplification of gravity waves is shown for a "potential barrier" ^[0(77)] = 
(a/a), for inflation (a = (rj /i]) ) followed by radiation and matter domination 3 (a = a rj(rj + fj ) ). 
At early times t < t\ , the scalar field is in its vacuum state (f)^ . At late times t w to, the scalar field 
will in general be described by a linear superposition of positive and negative frequency solutions of 
eq. (1) : 4>k(v) = a ^ ) i + ' ) + A mode with comoving wavenumber fc = (2na/\) is shown to leave 

the Hubble radius during inflation (t\), and re-enter it during matter domination (to)- (Figure not 
drawn to scale.) 

Fig. 2 

The ratio of the Quadrupole anisotropy from tensor and scalar waves is plotted against both 
spectral index n and p (in a oc t p ) for power law inflation. The maximum value p = 134 corresponds 
to n = 0.985 and is equivalent to ~ 67 e-foldings of slow-roll inflation. 

Fig. 3 

The biasing factor b\§ = (^)j" 6 , is shown plotted against the spectral index n, for three values 
of the Hubble parameter: H = 50, 75 & 100 kmsec" 1 Mpc~ x . The dotted lines correspond to b\% = 1 
and &i6 = 2.5. 

Fig. 4 

The COBE-normalised spectral energy density (in units of the critical density) and the amplitude, h x 
of gravity waves is shown (in figures a and b respectively) as a function of the wavelength ( and fre- 
quency /), for exponential inflation, quasi-exponential inflation (dotted line), and power law inflation 
aocf with p = 21 and p = 9. For comparison, the expected sensitivity of the Laser Interferometer 
Gravitywave Observatory (LIGO), (Christensen 1992) and of the projected "Beam in space" (space 
- interferometer) has also been plotted (Thorne 1988). 



